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Abstract. A relativistic radio component, which moves 
in a direction close to the sky plane, will increase in flux 
density when it decelerates. This effect is the basis for 
the qualitative model for GPS galaxies we present in this 
paper, which can explain their low-variability convex spec- 
trum, their compact double or compact symmetric mor- 
phology, and the lack of GPS quasars at similar rcdshifts. 
Components are expelled from the nucleus at relativistic 
speeds at a large angle to the line of sight, and are decel- 
erated (eg. by ram-pressure or entrainment of the external 
gas) before contributing to a mini-lobe. The young compo- 
nents arc Doppler boosted in the direction of motion but 
appear fainter for the observer. The non-relativistic mini- 
lobes dominate the structure and are responsible for the 
low variability in flux density and the convex radio spec- 
trum as well as the compact double angular morphology. 
Had the same source been orientated at a small angle to 
the line of sight, the young components would be boosted 
in the observer's direction resulting in a flat and variable 
radio spectrum at high frequencies. Hence the character- 
istic convex spectrum of a GPS source would not be seen. 
These sources at small angles to the line of sight are prob- 
ably identified with quasars, and are not recognized as 
GPS sources, but are embedded in the large population of 
flat spectrum variable quasars and BL Lac objects. This 
leads to a deficiency in GPS quasars at z < 1, but the 
model does not explain the population of GPS quasars at 
high rcdshift. However, there is increasing evidence that 
the high redshift GPS quasars are not related to the GPS 
galaxies. It is interesting for unification purposes to inves- 
tigate if deceleration is a general phenomenon in compact 
radio sources. We have therefore investigated if evidence 
for deceleration can be traced in flux density outbursts of 
highly variable radio sources, and have found that indeed 
the evolution of the spectral peak in such outbursts does 
suggest that deceleration may play an important role, at 
least in some of these sources. 



1. Introduction 

Gigahertz Peaked Spectrum (GPS) sources are a class of 
extragalactic radio source characterized by a convex radio 
spectrum peaking at a frequency of about 1 GHz. They are 
compact luminous radio sources. The turnovers in their 
radio spectra are believed to be due to synchrotron-self- 
absorption caused by the high density of relativistic elec- 
trons in these compact sources. The overall angular size 
of GPS sources is a function of their spectral peak: the 
lower the peak frequency the larger the radio source. GPS 
sources may therefore be related to the class of compact 
steep spectrum (CSS) sources, which have lower peak fre- 
quencies (~ 100 MHz) and larger angular and linear sizes. 

Optical identification programs carried out on GPS 
sources show the optical counterparts to be a mixture 
of galaxies and quasars (O'Dea 1991, Stanghellini et al. 
1993, de Vries 1995). The quasars tend to have very high 
redshifts, while the GPS-galaxies have redshifts which are 
generally well below unity (O'Dea et al. 1991). Within the 
context of the simple orientation unification scheme of ra- 
dio galaxies and quasars, it is puzzling that the redshift 
distribution of GPS galaxies is so different from that of 
GPS quasars. There is an increasing amount of evidence 
available in the literature that the GPS sources identified 
with galaxies and quasars are not related to each other and 
are different classes of objects; in addition to their differ- 
ent redshift distributions, it has been found that they have 
different distributions in rest-frame peak frequency, linear 
size and radio morphology (Stanghellini et al. 1996). 

Early VLBI data on GPS sources showed that many of 
them, especially those identified with galaxies, had a com- 
pact double (CD) morphology with components of simi- 
lar flux densities and spectral indices (Phillips and Mutcl 
1982). This led to the speculation that these are young 
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Fig. 1. 1144+352 at 5 GHz from the Caltech-Jodrell II survey 
(Henstock et al. 1995). The brighter component to the east has 
a flux density of 467 mjy and is mainly responsible for the ob- 
served flux variability over the last 20 years. The two brightest 
components separate from each other with an apparent veloc- 
ity of 1.2c/ft _1 (Giovannini et al. 1995). 



radio sources, in which the two components could be in- 
terpreted as being mini-lobes. More recent higher quality- 
images show that in many CDs weak central components 
with flatter spectral indices are found which are most nat- 
urally interpreted as being the central core. These sources 
were renamed compact symmetric objects (CSO), and it 
has been proposed that CSOs evolve into FRI and/or 
FRII sources (Readhead et al. 1994, Fanti et al. 1995). 
Stanghellini et al. (1996) has recently shown that about 
90% of the GPS sources identified with galaxies are CSOs. 

In this paper we present a model for GPS galaxies 
which can explain their spectral shape, their compact dou- 
ble or compact symmetric morphology, their apparent low 
flux density variability and the lack of GPS quasars at sim- 
ilar redshifts. In section 2.1, th e ob servational basis for the 
model is presented. In section |2.2[ the model is explained 
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and its implications for the appearance of GPS sources 
presented. In section ||, the possiblity that deceleration is 
a more general phenomenon is investigated, by studying 
the spectral evolution of flux density outbursts in highly 
variable radio sources. 



The Spectral and Angular Morphology of GPS 
Galaxies 



Fig. 2. 1946+708 observed at 5 GHz by Snellen et al. (in 
prep.). Taylor et al. (1997) found bi-directional motion for the 
inner components, while for the outer components no motion 
is detected. This indicates that the components are decelerated 
while moving outwards. 



2.1. The Starting Points 

The following three observations of properties must be 
explained by any model of GPS sources: 

1 The symmetric structures of GPS sources. The 

compact symmetric structure seen in most GPS galax- 
ies (Stanghellini et al, 1996) indicates that relativistic 
boosting is relatively unimportant, otherwise only ra- 
dio emission on the side directed towards us would be 
visible. This implies that these components are not 
ejected from the nucleus at small angles to the line of 
sight and at highly relativistic speeds. However it is 
not excluded that the components move at highly rel- 
ativistic speeds at large angle to the line of sight (see 
point 2). 

2 The relativistic velocities of components in 
nearby GPS galaxies. In two nearby GPS sources, 
1144+352 (Giovannini et al. 1995) and 1946+708 (Tay- 
lor et al. 1997), it has been shown that some compo- 
nents have relativistic motions. The source 1144+352, 
which is shown in figure |l|, consists of two compo- 
nents (ignoring faint large scale structure at long wave- 
lengths) which separate from each other with an ap- 
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parent velocity of \flcjhr 1 . The structure of 1946+708 
consists of a series of components (see figure ^ found 
by Taylor et al. to move away bi-directionally from a 
point between the two central components, probably 
the location of the nucleus. It is noteworthy that the 
components close to the centre move outwards faster 
than the components further away and that no mo- 
tion is found for the two components at the largest 
distances from the centre. The highest apparent veloc- 
ity measured in this source is 0.9c/ Ti" 1 (Taylor et al 
1997). 

3 The flux density variations of the nearby 
GPS source 1144+352. Although the GPS source 
1144+352 is variable on timescales of several years 
(Snellen et al. 1995), its peaked spectral shape does not 
change. The peak flux density of the source increased 
monotonically from 300 to 600 mJy over some 20 years, 
and is at the moment decreasing slowly again. In view 
of the simple structure of 1144+352, we can attribute 
this flux density variation completely to the eastern- 
most (brightest) component in this source (Giovannini 
et al. 1995), because it contributes 85% of the flux 
density at 5 GHz. 

The variability of 1144+352, the relativistic velocities 
of components in nearby GPS sources, and the large angles 
to the line of sight implied by symmetric morphologies 
lead us to the model discussed below. 

2.2. A Decelerated Component Model 

We shall model a GPS radio galaxy as a series of radio 
components expelled at different epochs, hence observed 
at different stages of their evolution. Each of those com- 
ponents is assumed to be expelled from the nucleus at 
relativistic speeds along a direction which makes a large 
angle to the line of sight. These components then decel- 
erate and expand as they interact with the surrounding 



medium. In section 3.1 we will go briefly into the possible 



nature of the expelled components and their deceleration 
mechanism, however this is not necessary for our qualita- 
tive approach and not important at the moment. 

Consider the evolution of one such component. Initially 
its velocity is highly relativistic (figure |3|, component 1) 
and the intensity of its radio emission is strongly Doppler 
boosted in the direction of motion. The component is as- 
sumed to be sufficiently compact to be optically thick at 
low frequencies as a result of synchrotron self absorption 
(SSA). Assuming that the component emits isotropically 
with a flux density Sq(v) at frequency i>, it will be ob- 
served to have a flux density of 



S(y&) = S (v)5 3 



(1) 



where 8 is the Doppler factor, 8 = (7(1 — (3 cos^)) -1 , 
j3 is the velocity in units of the velocity of light, 7 = 
1 / y/ 1 — (3 2 , and 9 the angle of the direction of the motion 



to the line of sight. For small values of 9, and highly rel- 
ativistic velocities this results in 8 >> 1. However, for a 
large angle to the line of sight, cos 9 will be small, hence 
8 w 7 _1 , which results in 8 < 1. As it evolves, the com- 
ponent will be decelerated to a modest relativistic speed 
(figure 0, component 2), as is seen in 1946+708, eg. by 
drag or entrainment of the external gas. The component 
also expands, resulting in a intrinsically lower SSA peak 
frequency and peak flux density. 

Figure ^ shows the Doppler factor as function of 7 for 
different angles to the line of sight. For small 9s and de- 
creasing (3, the Doppler factor decreases because 1 _p cosS 
decreases faster than I/7 increases. However at large an- 
gles to the line of sight, 8 increases with time. For 9 = 90°, 
the peak flux density, which is proportional to 8 3 , increases 
by a factor of 10 between (5 = 0.98 (7 = 5) and (3 = 0.90 
(7 = 2.3). Eventually the component is decelerated to 
barely relativistic speeds or slower and will contribute to 
a 'mini-lobe' (figure ||, component 3). 

We now consider the effect of the simultaneous pres- 
ence of the fast and slow components on the overall angu- 
lar and spectral morphology that would be observed in a 
population of radio sources (figure |J). 

i) For objects whose radio axes are oriented at a 
large angle to the line of sight, the slow, barely- or non- 
rclativistic components, which form mini-lobes on both 
sides of the nucleus, are dominant in both the radio spec- 
trum and angular morphology. The mini-lobes are still 
compact enough to be self absorbed, resulting in a 'giga- 
hertz peaked spectrum', and the radio source will have a 
compact symmetric morphology. The fact that 8 < 1 for 
the highly relativistic young components, makes them ap- 
pear faint. Therefore the flux density variability will be 
low due to the small contribution to the total flux density 
of the young and fast components; that is the creation of 
a new component will only have a small impact on the 
overall radio spectrum. 

ii) For objects whose radio axes are oriented at a small 
angle to the line of sight, the influence of Doppler boost- 
ing on the appearance of the (slow) mini-lobes is small. 
They still would be observed to have a 'gigahertz peaked' 
radio spectrum, were it not that the young components, 
which move towards the observer, are strongly boosted in 
flux density and blue-shifted in frequency and dominate 
the overall spectrum producing a peak at high frequency. 
The overall radio spectrum appears to be flat or inverted 
with a peak at high frequency (see figure ||) , and the ob- 
ject is not classified as a GPS source. Furthermore, the 
creation of new components produces a large effect on the 
radio spectrum, and results in the source being observed 
as (strongly) variable at high frequencies. 

According to the orientation unification scheme (eg. 
Barthel, 1989) a quasar whose radio axis is oriented at a 
large angle to the line of sight, is observed to be a radio 
galaxy. However in our model, a GPS galaxy observed 
at a small angle to the line of sight, would not result in 
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Fig. 3. The decelerated component model. The cross indicates the central engine. Components are expelled from the nucleus 
at relativistic speeds (1), then decelerate (2) and then contribute to a mini-lobe (3). The arrows indicate the velocity, and the 
grey ellipses the anisotropy of the observed radio emission due to Doppler boosting. 
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Fig. 4. The Doppler factor 8 to the third power, which is 
proportional to the peak flux density of a source, as function 
of 7 for an angle of 15°, 30° and 60° to the line of sight. 



a quasar with a GPS radio spectrum, but in a quasar or 
blazar with a flat or inverted variable spectrum. Therefore 
no population of GPS quasars at similar redshifts to the 
GPS galaxies will be observed. These quasars would hide 
as a subsample in the large population of flat spectrum 
quasars and BL Lac objects. 



3. The Evolution of a Single Component 

The physical properties of radio jets are widely discussed 
in the literature, but this has not yet converged to a single 
view. Fortunately, no detailed understanding of the un- 
derlying jet-physics is necessary to come to a qualitative 
description of the flux density evolution of a decelerat- 
ing radio component. However, we will discuss briefly the 
possible nature and deceleration mechanisms of the radio 
components in GPS sources. 

3.1. The Physical Nature of the Radio Components and 
the Deceleration Mechanism. 

In general it is assumed that radio jets can be treated 
as continuous fluids, but the jet-physics of GPS sources 
may be different from radio sources in general due to their 
young age; eg. a steady continuous jet may not have been 
formed yet. The most simple discrete radio source model 
is the adiabatic, spherically expanding source model (eg 
van der Laan 1966, Pauliny-Toth and Kellermann 1966), 
which was the first quantitative attempt to explain ra- 
dio outbursts. If the radio components could indeed be 
treated as individual 'plasmons' the deceleration mecha- 
nism can be easily explained due to ram-pressure of the 
surrounding medium. Such a scenario has been modeled 
by Christiansen et al. (1978). However, this scheme results 
in a large energy release by dissipation which should be 
observed in some way. 

In radio source models which consider jets to be con- 
tinuous, the radio components arc in general treated as 
shocks. In this case deceleration of the jets can be caused 
by entrainment of the external gas. When deceleration of 
jets was first discussed for FRI radio sources (eg. Begel- 
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Fig. 5. The decelerated component model combined with the orientation unification scheme. At a large angle to the line of 
sight, the Doppler boosting of component 1 does not make a significant contribution to the overall radio spectrum, resulting in a 
GPS spectral morphology. In this case the object is seen as a galaxy. At a small angle to the line of sight, the young components 
moving towards us are strongly Doppler boosted, resulting in a one sided jet morphology and a variable flat spectrum. In this 
case the object is seen as a quasar. Hence the lack of GPS quasars at similar redshifts to the GPS galaxies is due to their not 
being recognized as GPS sources, since they are embedded in the large population of variable quasars and BL Lac objects, with 
radio emission only at sub-arcsecond scale. 



man 1982, Scheuer 1983), it was pointed out that if an ini- 
tially relativistic jet decelerates to sub-relativistic speeds 
much of its kinetic energy is thermalized, but that the jet 
can be recollimated by an external pressure gradient. More 
detailed analyses have been carried out by Bicknell (1994, 
1995), Komissarov (1994). Bowman et al. (1996) claim 
that the energy released by dissipation is re-converted 
continuously into kinetic energy since the jet is propa- 
gating down a pressure gradient, and that deceleration 
without catastrophic kinetic energy loss is relatively easy 
to achieve for relativistic jets compared to non-relativistic 
ones. Recently, de Young (1997) concluded that for a con- 
stant hot-spot advance speed for CSOs as derived from 
VLBI data by Readhead et al. (1996) it is necessary to 
have entrainment along most of the jet, therefore allowing 
deceleration of the jet. 



3.2. The Qualitative Flux Density Evolution of a Deceler- 
ating Component 

If the deceleration model described above is valid, then 
deceleration should have a large influence on the spectral 
evolution of flux density outbursts observed in quasars and 
BL Lac objects. It is therefore interesting to investigate 
these flux density outbursts, and see if there is any sig- 
nature of deceleration. First of all we have to discuss how 
flux density outbursts are explained in current models. 

The adiabatic expanding source model agreed well 
with the observed spectral evolution at low frequencies, 
but it was less successful in explaining the high-frequency 
behaviour due to the earliest stages of the outbursts. 
Closer agreement was obtained using shock models, in 
which the radio outbursts are explained as due to prop- 
agation of shocks in relativistic jets (Marscher and Gear 
1985, Hughes et al. 1989). The general applicability of the 
shock models is now widely accepted, although there are 
some features that these models cannot yet explain prop- 
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Fig. 6. The evolution of the SSA spectral peak of a radio 
outburst according to the shock model (Marscher and Gear 
1985). Note that in the first phase, where inverse Compton 
losses are dominant, the peak frequency is expected to decrease 
with time. 



Fig. 7. The evolution of the SSA spectral peak of a radio 
outburst according to the decelerated component model. Note 
that in the first phase, where deceleration is dominant, the 
peak frequency is expected to increase with time. 



erly (eg. the growth stage of the shock). Three phases 
can be distinguished in the shock model (see figure ||). fn 
the first phase, Compton losses are dominant, resulting in 
a dependence of the SSA peak flux density S pea k to the 
peak frequency v pea k of the shocked component given by 
Speak oc v p ^[ k 2 . In the second phase, synchrotron losses 
are assumed to be dominant, resulting in S pea k oc v~® a ®° ■ 
In the last phase, expansion losses are dominant, which 
is comparable to the adiabatic expanding radio source 
model, resulting in S pea k oc v peak - Note that the exact 
relations between S pea k and v pea k in this model are both 
dependent on the relativistic electron energy distribution 
and the change in magnetic field strength along the jet (see 
Marscher and Gear 1985). The predictions of this shock 
model agree well with spectral observations (eg. Lainela 
1994). 

What influence can deceleration have on the spectral 
evolution of a radio outburst? Assuming that deceleration 
is the dominant factor, the change in the Doppler factor 6 
results in a change of the SSA spectral peak. From equa- 
tion 1 we find that AS pea k oc (A(S) 3 and /S.v pea k oc S, and 
therefore that AS pea k oc (Av pea k) 3 . If the velocity of the 
component is highly relativistic, and 7 >> 1, deceleration 
will cause an increase in <5, resulting in an increase in peak 
flux density and peak frequency, even for a small angle to 
the line of sight (see fig. |). For example, at an angle of 
20° to the line of sight, 6 is maximized at 7=2.9. Hence for 



7 > 2.9 deceleration will cause an increase in 8. Between 
7 = 10 and 7 = 2.9, the observed peak flux density will in- 
crease by a factor of 40, and the peak frequency increases 
by a factor of 3.4. During the deceleration of a compo- 
nent, it will also expand. At a certain point, the change in 
S will be small compared to the changes in the spectrum 
due to the expansion. This period will be similar to the last 
phase in the general shock model. The effect of combined 
deceleration and adiabatic expansion is shown in figure 7, 
in which it is clear that the predicted spectral evolution 
is comparable to that expected for the generalized shock 
model, except in the initial phase of the outburst. The de- 
celeration phase in fact replaces the inverse Compton and 
synchrotron phases. How can we distinguish between the 
two models? A detection of an increase in peak frequency 
combined with an increase in peak flux density of a flux 
density outburst would underline the importance of decel- 
eration. However, this is difficult to observe, because it is 
expected to occur in the very early phase of the outburst. 
We have searched the literature for examples, as we now 
describe. 



3.3. Observations of the Spectral Evolution of Flux Den- 
sity Outbursts. 

To be amenable for investigation of its spectral evolution, 
a radio outburst should have properties which obey the 
following criteria: 
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1) The outburst should be a single event uncontaminated 
by previous outbursts. Blending of spectral features 
due to different outbursts make detailed investigation 
impossible. 

2) The radio spectrum during the initial period of the out- 
burst, before maximum flux density is reached, must 
be well sampled. This implies that the timescale of the 
outburst should be sufficiently long for this to be ac- 
complished. 

3) The radio source must be monitored at three or more 
frequencies to obtain sufficient information about the 
evolution of the spectral peak. Furthermore the highest 
frequency at which the source is monitored should be 
in the optically thin part of the spectrum, to allow 
the peak to be properly fitted. This implies that the 
maximum flux density of the outburst at the highest 
observed frequency should be lower than the maximum 
flux density at lower frequencies. 

For many years large dedicated monitoring programs 
to study the time variability of radio sources have been 
carried out by the University of Michigan at 4.8, 8.0 and 
14.5 GHz (Aller ct al. 1985) and the Metsahovi Observa- 
tory at 22 and 37 GHz (Saloncn 1983, 1987, Terasranta 
1992). We used data obtained in the course of these 
two monitoring programs to search for radio outbursts 
with properties obeying the criteria mentioned above. 
The best cases we found are the outbursts in 3C345 and 
3C454.3 peaking in 1981 and 1982 respectively. Further- 
more, Margo Aller provided us with unpublished data 
from an outburst in 3C454.3 in 1986. All three outbursts 
have well sampled observations at the 5 frequencies men- 
tioned above. The contribution of emission not related to 
the outburst, often refered to as quiescent emission, was 
estimated from the flux densities in the period before the 
outburst, and then subtracted from the total observed flux 
densities. For 3C345 the quiescent level is about 7.0 Jy at 
all observed frequencies, and for 3C454.3 the quiescent 
emission was estimated to be S q = 10. Of -0 - 2 . The result- 
ing spectra were fitted with a least squares fitting routine, 
in which the optically thin spectral index was fixed at 
a = —0.75. Note that if the optically thin spectral in- 
dex were to change during the outburst, it would become 
steeper in time (due to aging of the high energy electrons) 
and results in an estimation of the peak frequency which is 
too low. Therefore keeping a constant cannot be the rea- 
son for finding an increase in peak frequency with time. 
The resulting evolution of the spectral peak in time for 
the three different outbursts are shown in figure 7. In the 
outburst in 3C345 and the second outburst in 3C454.3, 
the peak frequency does not significantly change in the 
initial stage of the outburst. However the first outburst 
of 3C454.3 seems to have an increase in peak frequency 
in the initial stage of the outburst, which is expected if 
deceleration plays an important role. In all three cases 
the peak frequency to peak flux density relation is steeper 



than expected for the general shock model (S pea k ~ v pea k^ 
with k = —5/2), we find k = —6.7 for the 3C345 outburst 
and k = +1.9 and +3.6 for the 3C454.3 outbursts. Hence 
these observations support the idea that deceleration plays 
a role in the spectral evolution of radio outbursts and, in 
some cases, may have a dominant role, e.g. 3C454.3. 

Another method of investigating the spectral evolution 
of a radio source during an outburst would be to moni- 
tor its spectral index between two frequencies far enough 
apart that the spectral peak is enclosed between them. 
This has been done by Stevens et al. (1995) for the blazar 
PKS 0420-014. They have monitored the spectral index 
between 90 and 22 GHz and find clear relations between 
the 90 GHz flux density and the spectral index for the ris- 
ing and falling phases of its radio flux density evolution. 
These relations can be converted to a relation between 
peak flux density and peak frequency for both stages, as- 
suming the quiescent level (given by Stevens et al.) and a 
homogeneous self absorbed source with an optically thin 
spectral index, for which we take the value of Stevens et 
al. of -0.44 and assume it to be constant. During the rise 
in flux density the evolution of the spectral peak obeys 
the relation S pea k oc Vp eak , where k = +2.0to.7- During 
the decay in flux density it follows the same relation with 
k = +1.1 + 0.1. Hence these observations also support the 
idea of deceleration. 

4. Discussion 

The decelerated component model which we discuss in this 
paper can explain several properties of Gigahertz Peaked 
Spectrum radio galaxies. While consistent with the typi- 
cal two-sided morphology and low variability, it can also 
explain the lack of GPS quasars at similar redshifts as the 
GPS galaxies, and account for the spectral evolution of 
flux density outbursts in highly variable radio sources. In 
this scheme, 1946+708, our prototype compact symmet- 
ric GPS source, is oriented close to the plane of the sky 
and would not have had a gigahertz peaked spectrum, if it 
were oriented at a small angle to the line of sight. It has a 
peak frequency of about 2 GHz and a peak flux density of 
about 1 Jy, and is believed to lie at a large angle to the line 
of sight, 9 = 75° (Taylor et al. 1995). The measured flux 
density of the innermost component is about 50 mJy at 5 
GHz. If the angle to the line of sight had been 30°, and 
assuming that its velocity is 0.9c, the doppler factor had 
been a factor 3.5 higher and this young component would 
contribute 2 Jy at 17.5 GHz. Hence the overall spectrum 
of this radio source would appear flat above 1 GHz and 
it would not be classified as a GPS source. Furthermore 
the rapid change in Doppler boosting and expansion of 
the young component would make the radio source ap- 
pear significantly more variable. For objects whose radio 
axes are oriented at small angle to the line of sight, the 
young components will be dominant at high frequencies, 
but at about 1 GHz, they would be observed in the opti- 
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Fig. 8. The evolution of the spectral peak during outbursts in 3C345 (left) and 3C454.3 (middle and right). The numbers 
indicate the dates of the observations. The dates for 3C345 are 1980.5 and every half a year after. The dates for 3C454.3 are 
1981.1 and every 0.3 years later for the first outburst, and 1985.5 and every half year later for the second outburst. Note that, 
especially the 1982 outburst of 3C454.3, undergoes an increase in peak frequency during the first stage of the flare which is not 
expected for the generalized shock model. However this is expected if deceleration is important. 



cally thick part of their spectra where Doppler boosting is 
only mildly important (S oc <5 5 ). Hence their 1 GHz flux 
densities would be similar to those observed if they were 
oriented at a large angle to the line of sight. A GPS galaxy 
observed at 9 < 20° will appear as a flat spectrum vari- 
able quasar or a Bl Lac object. In a flux density limited 
sample at low frequency, only 1 flat spectrum quasar with 
no extended emission is expected per 20 GPS galaxies as- 
suming a parent population with a random distribution of 
0. 

This corresponds to only 1 to 2 objects with a flux 
density greater than 1 Jy at 1 GHz. An example of such 
an object could be PKS 1413+135 (Pcrlman et al. 1994), 
which has no measured extended emission, a "flat" vari- 
able spectrum at high frequency, and at VLBI scales a 
component opposite to the jet which can be interpreted 
as a non-Dopplcr-boosted mini-lobe on the receding side 
of the source. Another candidate GPS source seen at a 
small angle to the line of sight is 1504+377 (CJI Survey, 
Polatidis et al. 1995). 

Although the angle to the line of sight has to be large 
to produce a Gigahertz Peaked Spectrum, will not often 
be exactly 90°. Therefore one side will have a component 
directed towards us and the other side will have a com- 
ponent directed away from us. This will result in different 
Doppler boosting for the approaching and receding jet, 
especially for the inner part of the jet which is moving 
the fastest. Such a situation could account for the mor- 
phologies observed in sources like 0108+388, 0404+768, 
2352+495 and 0710+439 (Taylor et al. 1996): While the 
outer components are comparable in flux density, the in- 
ner components are only visible in the approaching jet. 
An exception to this is our prime example 1946+708. The 



inner two components are comparable in flux density for 
both the approaching and receding jet. This could be due 
to the fact that we observe these components in the op- 
tical thick, inverted parts of their spectra which are not 
very sensitive to Doppler boosting. However the outer two 
components, which show no motion, differ by a factor of 
~ 5 in flux density. This could be caused by inhomogencity 
of the medium for example, which is proposed by Taylor 
et al. 1995. However, there is a possibility that due to 
the fact that the approaching component is observed at 
an intrinsically later epoch than the receding component, 
that this difference in flux density is a result of observ- 
ing the components at different evolutionary stages. VLBI 
observations at different frequencies are needed to further 
investigate this. 

So far, we have ignored the Gigahertz Peaked Spec- 
trum radio sources identified with quasars, which are 
found at high redshift. The model does not explain their 
existence. However, the GPS radio galaxies and quasars 
have a different distribution in redshift, rest frame peak 
frequency, linear size and VLBI morphology, and are prob- 
ably different classes of objects (Stanghellini et al. 1996). 
Still it is confusing that these quasars, which are believed 
to be observed at a small angle to the line of sight, do not 
undergo significant Doppler boosting. 

We have shown that deceleration could play an impor- 
tant role in 3C454.3 and 0420-014 and have hypothesized 
that it plays an important role in GPS galaxies and in 
other flat spectrum variable quasars and BL-Lacs. We also 
have noted that there is in fact very little data available in 
which the spectral evolution of individual outbursts can be 
unambiguously followed to look for evidence of decelera- 
tion. However it remains true that these outbursts can be 



Snellen et al.: Deceleration and the Morphologies of GPS sources 



9 



= 5° 



6 = 2.5° 



d 
o 



5-i 
O 
'A 
O 



5m 



10 




1 



10 



Doppler Factor 5 



Fig. 10. The influence of 7 and the angle to the line of sight 8 on the Doppler factor S and the apparent proper motion f3 apP - 
The solid lines represent a range of 7s for a certain fixed 8, while the dotted lines represent a range of 8s for a certain 7. The 
areas A and B indicate the parameter space for which an increase in S can be due to deceleration and acceleration respectively. 



exceptions. Furthermore, the quiesent radio flux density 
was assumed to be constant, and the change in radio flux 
density is all contributed by one component. This does 
not have to be the case, and the influence of other compo- 
nents can change the observed relation between peak flux 
density and peak frequency. 

Ideally, one should analyse an outburst with VLBI at 
multiple frequencies. A first attempt has been carried out 
by Lobanov and Zensus (1996) using data on 3C345 at 8 
different frequencies ranging from 1.4 to 89 GHz. Indeed a 
clear increase in peak frequency for two components was 
found which is in contradiction with the basic shock mod- 
els. Lobanov and Zensus also conclude that an increase in 
the Doppler factor S is the probable cause; however they 
believe this is due to acceleration and/or a change in the 
angle to the line of sight 9. There is a hint for acceleration 
in the measured angular velocities for these components 



(Zensus, Cohen and Unwin 1995), but the evidence for 
this is not convincing. If a change in the Doppler factor d 
causes the increase in peak frequency (and peak flux den- 
sity) in 3C345, it does not have to result in a change in 
the observed proper motion. In figure [To], the influence of 
7 and the angle to the line of sight 8 on both the Doppler 
factor S and the apparent proper motion f3 app is shown. 7 
varies for the solid lines while 9 varies for the dotted lines. 
The grey areas A and B indicate the regions in parameter 
space for which deceleration and acceleration respectively 
can cause the increase in Doppler boosting. If you follow 
a solid line in region A, S increases while 7 decreases, the 
reverse is true in region B. In region A the change in ap- 
parent proper motion is only small, while in region B, 
Papp increases significantly. Moreover in region A, the an- 
gle to the line of sight 9 is in the range 10° — 45°, while 
for region B, 9 is in the range 0° — 10°. Hence, due to 
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Fig. 9. The evolution of the peak frequency and peak flux 
density of a flux density outburst of 0420-014, derived from 
the spectral index monitoring by Stevens et al. (1995) at 90 
and 22 GHz, assuming a homogeneous self-absorbed source. 
During the rise in flux density the evolution of the spectral 
peak the peak frequency increases, hence these observations 
also support the idea of deceleration. 

retardation effects, an outburst undergoing an increase in 
S due to acceleration should be observed to have a much 
smaller timescale than an outburst undergoing an increase 
in S due to deceleration. The sources discussed in this pa- 
per are selected for their slow and well studied variability, 
and are therefore most likely not at the smallest angles to 
the line of sight compared to other highly variable radio 
sources. If this decelerated component model is correct, 
then no increase in peak frequency is expected for the 
rapidly variable sources with 9 < 10°, unless 7 > 20. Note 
that a change in 9, if it is combined with the right 7, 
can also cause S to increase without changing (3 app signif- 
icantly. However this is especially the case for small 9, for 
which only a small change in 9 can already significantly 
increase S. 

5. Conclusions 

A qualitative model for GPS galaxies has been presented, 
in which components are expelled from a nucleus at rcl- 
ativistic speeds, decelerate and possibly contribute to a 
mini-lobe. At large angles to the line of sight, the minilobes 
are dominant because the Dopplcr factor of the young 
components is smaller than unity along the line of sight 
direction. At small angles to the line of sight these young 



components appear much brighter and become dominant. 
This model is consistent with the symmetric angular 
morphology of the mini-lobes of GPS galaxies and their 
low variability. It can explain the apparent lack of GPS 
quasars at similar redshifts to the GPS galaxies. Further- 
more it provides an alternative explanation for the evolu- 
tion of the spectral peak of flux density outbursts in highly 
variable radio sources. 
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